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Abstract 



Crossover motifs are integral components for designing DNA based nanos- 
tructures and nanomechanical devices due to their enhanced rigidity com- 
pared to the normal B-DNA. Although the structural rigidity of the double 
helix B-DNA has been investigated extensively using both experimental and 
theoretical tools, to date there is no quantitative information about struc- 
tural rigidity and the mechanical strength of parallel crossover DNA motifs. 
We have used fully atomistic molecular dynamics simulations in explicit 
solvent to get the force-extension curve of parallel DNA nanostructures to 
characterize their mechanical rigidity. In the presence of mono-valent Na"*" 
ions, we find that the stretch modulus (71) of the paranemic crossover (PX) 
and its topo-isomer JX DNA structure is significantly higher (~ 30 %) com- 
pared to normal B-DNA of the same sequence and length. However, this 
is in contrast to the original expectation that these motifs are almost twice 
rigid compared to the double-stranded B-DNA. When the DNA motif is 
surrounded by a solvent with Mg^"*" counterions, we find an enhanced rigid- 
ity compared to Na"^ environment due to the electrostatic screening effects 
arising from the divalent nature of Mg^+ ions. This is the first direct de- 
termination of the mechanical strength of these crossover motifs which can 
be useful for the design of suitable DNA for DNA based nanostructures and 
nanomechanical devices with improved structural rigidity. 

Key words: PX/JX DNA motif; Atomistic Simulations; Force-extension 
Curves; Stretch Modulus; Torsion Angles 
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1 Introduction 

The structural properties that enable DNA to serve so effectively as genetic 
material can also be used for DNA based computation and for DNA nan- 
otechnology ([l|). The success of these applications to large extent, depends 
on the complementarity that leads to the pairing of the strands of the DNA 
double helix. This pairing can be exploited to assemble more complex mo- 
tifs, based on branched structures with sticky ends that are promising for 
building macromolecular structures ([l|, 0) . The combination of sticky ended 
ligation i) and stable branched BNA species (jj) has permitted investiga- 
tors to construct DNA molecules whose edges have connectivity of a cube (jSl) 
and of a truncated octahedron (0) . Similarly, branched DNA molecules have 
provided the basis for the syntheses of variety of structures such as knots 
|3, @) , Borromean rings (0) and DNA tweezers (fiol ) . One of the early goals 
of DNA nanotechnology was to construct precisely configured materials on 
a much larger scale ([l|, llll ). However the key feature that was lacking in 
such constructions were rigid DNA molecules, as flexible components failed 
to maintain the same spatial relationships between each member of a set 
(fl^ ). DNA anti-parallel double crossover molecules (DX) (0), analogous 
to intermediates in genetic recombination, were found to be about twice as 
rigid as a normal DNA molecules (114^ a nd were subsequently used for cre- 
ating two dimensional DNA lattices (|l5l ). 



The concept of double crossover molecules stems from the concept of 
formation of HoUiday junction (0), which is the most prominent interme- 
diate in genetic recombination. These crossover molecules usually contain 
two Holliday junctions connected by two double helical arms. Depending 
on whether the crossovers are between strands of the same (called parallel) 
or opposite (called anti-parallel) polarity, as well as the number of dou- 
ble helical half-turns (even or odd) between the two crossovers, different 
double crossover molecules are generated. It is known that anti-parallel 
double helical molecules are better behaved than the parallel ones and as 
mentioned above these anti-parallel molecules have found numerous appli- 
cations in nano construction The concept of crossover molecules also 
led to the discovery of paranemic crossover molecules (PX) and its topoi- 
somers, the various JX molecules. The PX motif arises from the fusion of 
two parallel double helices by reciprocal exchange at every possible point 
where the two helices come in contact, whereas its topoisomer JX, contains 
i adjacent sites where backbones juxtapose without crossing over (fisl ). The 
rigidity of these crossover motifs has led to wide applicability in DNA based 
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nanotechnology (fl9l ). 



Nevertheless, there is lack of quantitative knowledge about the rigidity of 
the paranemic PX DNA and its topoisomers JX DNA, even though they find 
applications in nano constructions, owing to their inherent rigidity. Using 
atomistic simulation studies on PX and JX structures we have shown that 
among all synthesized PX and JX DNA structures, PX-6:5 motif is the most 
stable and in comparison with JX motifs the rigidity order runs as PX-6:5 



> JXi > JX2 > JX3 > JX4 ([201, 121I ) . Our predictions on thermodynamic 



calculations are in consistence with the recent experimental studies using 



calorimetric and denaturing gradient gel electrophoretic methods (l22l ). In 
this report we make an attempt to get a quantitative measure of the rigidity 
of these structures from their force-extension curves. The concept that has 
been adapted is to pull these structures with an external force in constant- 
force ensemble (0) and from the internal energy as a function of extension 
of the DNA motif, we obtain the stretch modulus using Hooke's law of elas- 
ticity. The stretch modulus in turn gives an estimate of the rigidity of the 
motifs. 

Experimental and theoretical studies have established that when sub- 
jected to an external force, double-stranded DNA (dsDNA) exhibits different 
force-extension regimes, stemming from its unique double- helix structures 



(|24l427l). For example, in the low force regime, the elasticity of dsDNA is 
entropy dominated and the experimental force-extension data obtained with 
the magnetic bead method can be described very well by the standard en- 
tropic worm-like chain model (24, 28, 2^); in the high force regime (starting 
from 10 pN to 70 pN), which is accessible with optical tweezers or atomic 
force microscopy, when the external force is comparable with the basepair 
stacking interaction in dsDNA the polymers can be suddenly stretched to 



about 1.7 times its B-form length in a very narrow force range (|30l . l3ll). An 
explanation of this regime is attributed to the short-range nature of base- 
pair stacking interactions. Typically in this regime, the increase in length of 
the DNA molecule is due to large distortion in the double helical structure 



resulting in a ladder-like structure which is termed as S-DNA (j3ll 434l) where 
most of the H-bonds involving basepairing are intact. It has been shown that 
depending on the experimental conditions such as the force attachment and 
salt concentration, dsDNA can undergo either strand-unpeeling transition 
to ssDNA a a, or transition to S-DNA where basep^irs are intact (or 



B-to-S transition) ([3l|, 1371 ). In the elastic regime of force-extension curve. 



DNA behaves similar to that of a spring and follows the simple Hooke's Law: 
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F = YAAL/L where Y is the Young Modulus and A is the cross-sectional 
area of the DNA. The product of Young modulus and area, Y x A gives 
stretch modulus, 71. At large forces, the stacking potential can no longer 
stabilize the B-form conformation of dsDNA and the (optimally) stacked 
helical pattern is severely distorted (jssl ). and therefore a structural transi- 
tion from canonical B-form to a new overstretched conformation. While all 
these features are well established for dsDNA and ssDNA, no theoretical 
or experimental studies exist for the force-extension behavior of cross-over 
DNA molecules such as DX, PX and various JX structures. This paper is 
the first attempt to address such lacunas. The paper is organized as follows: 
in the next section we describe building of PX/JX nanostructures and the 
simulation methodologies. In section [3] we describe the results from our fully 
atomistic DNA pulling simulation and in section |4] we give a summary of 
the force-extension results and conclude. 



2 Computational Details and Methodology 
2.1 Building PX/JX structures 

In order to characterize the stiffness of various cross-over DNA motifs, the- 
oretically we have tried to mimic nano-manipulation techniques like AFM, 
magnetic tweezers, or laser tweezers in our fully atomistic simulations. The 
basic protocol in our simulation technique was to put the DNA structure in 
counterion solution, keep one end of the structure anchored, and pull the 
other end along the PX molecule dyad axis (parallel to the helix axes) with 
an external force and determine the force-extension curves. The various 
DNA structures studied are PX-6:5, JXi, JX2, AT-6:5, GC-6:5 and normal 
B-DNA which are four-strand complexes of DNA paired with one another. 
Initial structures of PX-6:5, JXi and JX2 were shown in Fig. [H We have 
also studied the force-extension behavior of a double strand DNA having 
the same length (38 bps) and the sequence of one of the helices of the PX 
structure. This will allow us to compare the force-extension behavior of the 
normal double strand B-DNA with that of the cross-over DNA nanostruc- 
tures. The procedure for constructing these structures is as follows: 

a. Building the DNA double helices: Regular B-DNA has between 10 and 
10.5 base pairs per helical turn. Hence by varying the twist angle of 
a selected number of base pairs, we can create B-DNA structures with 
between 8 and 18 base pairs per helical turn. Table 1 of (0) shows the 
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twist angles used for building the various PX structures. We assigned 
the same twist angle for all the base pairs in the helical half turn. The 
helical rise value of 3.4 A was used to build the PX structures. 

b. Building the crossover points: When a double helix is built in Namot2 
(version 2.2.) (jsol ). the molecules are oriented so that the 5' and 3' 
ends of the double helices are parallel to the y-axis. To create realistic 
crossover structures, it is necessary to rotate the individual helices so 
that the desired crossover points are closest to each other (rotation an- 
gles shown in Table 2 of (|2d )). To find this point we wrote a program 
that starts with the first crossover point and rotates the first helix in 
1° increment to find the rotation leading to the shortest distance be- 
tween these crossover points. Once found, the first helix is rotated by 
the prescribed value and held steady while the second helix is rotated 
and the shortest distance between the crossover points determined. The 
second helix is rotated 180° more than the 1st helix so that the helices 
are arranged as shown in Fig [TJ The crossovers were then created using 
the "nick" and "link" commands in Namot2. These structures are saved 
in the Protein Database (PDB) file format. Fig. [1] shows the snapshot 
of the built PX, JXi and JX2 structures after minimizations. 



2.2 Simulation Details for the PX/JX structures 



All molecular dynamics simulations were carried out using AMBERS suite 



of programs MQ') using the Amber 2003 (ff03) force fields (l4ll. |42 ) a nd the 



4J) inter- 



TIP3P model (j43l ) for water. For Mg^"*" we have used the Aqvist 
action parameters. This initial crossover DNA motif is then solvated with 
TIP3P model (0) water box using the LEaP module in AMBERS (0). 
In addition, some water molecules were replaced by Na"*" (Mg^"*") counte- 
rions to neutralize the negative charge on the phosphate backbone groups 
of the PX/JX DNA structure. We have used Aqvist parameter f|44) set 
to describe the ion-water as well as ion-DNA interactions. Aqvist (j44l ) pa- 
rameters reproduce the location of the first peak of the radial distribution 
function for the corresponding ion-water and the solvation free energies of 
several ionic species. After neutralizing the system with counterions, the 
concentration of Na+ (Mg2+) in the PX-6:5, JXi, JX2, AT-6:5, GC:6-5 and 
normal B-DNA crossover systems is 246 mM (160 mM) where as for a double 
helix B-DNA, the concentration of Na+ (Mg^+) is 254 mM (140 mM). The 
ion specificity and electrostatic interactions play crucial in the DNA func- 
tioning and protein-DNA binding mechanism (j45l . l46l ). The LEaP module 
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works by constructing Coulombic potential on a grid of 1 A resolution and 
then placing ions one at a time at the highest electrostatic potential. Once 
the placements of all the ions are done using the above method long MD 
simulations ensure that they sample all the available space around DNA 
and preferentially visit electronegative sites. In fact the initial placement 
of ions should not influence the final results provided long simulation was 
performed before the pulling runs were done. The radial distribution func- 
tion provided in the Fig. SI and S2 of supplementary material also reveal 
this fact that the counterions remain associated in the close proximity (~ 
10 A) of the DNA. We have also made sure that no counterion is "stuck" 
to the DNA molecule anywhere and observed diffusive behavior. Hence- 
forth we believe that the initial ion distribution is important to study DNA 
properties but ultimately equilibration is the key for avoiding such initial 
ion position dependency. The box dimensions were chosen in order to en- 
sure a 10 A solvation shell around the DNA structure in its fully extended 
form when the DNA is in overstretched regime. This procedure resulted 
in solvated structures, containing 99734 atoms for PX-6:5, JXi and JX2; 
97326 atoms for 38-mer B-DNA; 99781 atoms for PX-6:5 structure having 
only AT sequence and 99705 atoms for PX-6:5 structure having only GC 
sequence in a box of dimensions 54 x 77 x 299 A^ when Na"*" ions are present 
as counterions. Translational center of mass motions were removed every 
1000 MD steps. The trajectory was saved at a frequency of 2 ps. We have 
used periodic boundary conditions in all three directions for the water box 
during the simulation. Bond lengths involving bonds to hydrogen atoms 



were constrained using SHAKE algorithm (j47l ). During the minimization, 
the PX/JX structure was fixed in their starting conformations using har- 
monic potential with a force constant of 500 kcal/mol-A^. This allowed the 
water molecules to reorganize to eliminate bad contacts within the system. 
The minimized structures were then subjected to 40 ps of MD, using 1 fs 
time step for integration. During the MD, the system was gradually heated 
from to 300 K using weak 20 kcal/mol-A^ harmonic constraints on the 
solute to its starting structure. NPT-MD was used to get the correct solvent 
density corresponding to experimental value of the density. Lastly, we have 
carried out pulling the motif in NVT MD with 2 fs integration time step 
using a heat bath coupling time constant of 1 ps. 

The external pulling force was applied at one end on 03' and 05' atoms 
on each strand and the other 05' and 03' atoms of the strands on the other 
end were fixed with large force constant of 5000 kcal/mol-A^. The force 
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was applied along the direction of the end-to-end vector joining 03' and 05' 
atoms. The external force started at pN and increased linearly with time 
up to 1000 pN with a rate of force of 10~^ pN/fs. For comparison, we have 
also done the simulation with faster pulling rate at 10~^ pN/fs. It should 
be mentioned that the typical pulling rate in an AFM experiment is of the 
order of 10^ pN/s whereas the slowest pulling rate achieved in our simula- 
tion is 10^^ pN/s (or 10~^ pN/fs) due to computing limitations. Theoretical 
model suggests that with the increased rate of force, the dynamic strength 
of the molecular adhesion bonds increases (0). Therefore, since our sim- 
ulation employs higher rates of force compared to AFM pulling rates, the 
calculated stretch modulus will be higher than that of the stretch modulus 
calculated in an AFM pulling experiment. Our pulling protocol was vali- 
dated for B-DNA (1231 ) and is expected to be applicable for PX/JX motifs 
also. In Fig. S3(a) we have shown the force-extension curve that is consis- 
tent with experimental curve. Note that the plateau region was observed at 
95 pN instead of 65 pN which has its origin in fast rate of pulling which is 
inherent to computer simulations. This MD simulation has been done with 
a pulling rate of 10"^ pN/fs, the slowest pulling rate we could achieve as 
of now. We have also done the pulling simulation at various rates. In Fig. 
S3(b) we plot the plateau force as a function of the pulling rate. From this 
plot we can see that by extrapolating to slower pulling rates, one can get 
a plateau region like that of AFM pulling experiments. Hence we expect 
that we can exactly get the experimental force-extension curve when pulled 
much slower than 10~^ pN/fs such as in AFM or optical tweezers. Apart 
from the rate dependency, the stretch modulus calculation is in accordance 
with the experimental result. 



3 Results and Discussion 

3.1 Force-extension behavior 

Fig. [2] gives the force-extension curve for PX-6:5 DNA motif in the pres- 
ence of Na"*" counterions. The force-extension curve consists of an entropic 
regime where, the extension of DNA beyond its contour length is negligi- 
ble and this regime continues until 150 pN. Beyond this forces and up to 
200 pN, DNA extends to about 30 % with most of the H-bonds are still 
intact. It has been shown that depending on the experimental conditions 
such as the force attachment and salt concentration, dsDNA can undergo ei- 



ther strand-unpeeling transition to ssDNA (35l. l36l). or transition to S-DNA 
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where basepairs are intact (or B-to-S transition) (j3ll . l37l). This is foUowed 
by an overstretching plateau region resulting in an extension of 1.7 times 
the initial contour length where the DNA helical structure starts to deform. 
Beyond this plateau regime is the overstretched structure of DNA which 
is followed by strand separation. It is worth mentioning that at low rates 
of force, the above mentioned regimes in the force-extension curves shifts 
to much lower force values. For comparison we have also shown the force- 
extension behavior of the normal double strand B-DNA which has same 
length and sequence as one of the helices of the PX-6:5 structure. From the 
force-extension curve we can calculate the stretch modulus of DNA in the 
elastic regime using Hooke's Law. However, it is known from previous theo- 
retical and experimental studies on the DNA force-extension behavior that 
at high force regimes a conventional WLC model and its variant does not 
reproduce the force-extension behavior very well. Previous theoretical stud- 
ies have demonstrated that at high force regime, non-linear elasticity plays 
an important role governing the force-extension behavior of the DNA (|49l ) 
and the WLC model is inadequate to describe the force-extension behavior 
(0, ISlI ). Therefore, we have also used the following polynomial function to 
fit the simulated force-extension curve for whole force regime 

n — 1 



n=l ^ ^' ^ n=l 



Integrating the above equation, we get energy 

oo / r \ " °° 



n \Ln I ^-^ n 

n=-z ^ " ^ n=2 

Where £"0 is the energy of the DNA when the extension is zero, L is the 
length of the DNA and Lq is the initial contour length of the DNA. The coef- 
ficients 7n for various values of n gives various elastic modulii. For example, 
71 gives the linear stretch modulus. Stretch modulus 71 can be calculated 
from both the above equations by fitting the force and energy as a function 
of the applied strain. From the above two equations, it can be seen that 
the linear term coefficient and the quadratic term coefficient will give the 
stretch modulus, respectively. We have done a sequential fit to the energy 
vs strain curve similar to ref (fiol ) where the fitting is done for the leading 
quadratic form of energy vs strain plot for restricted data that are taken 
around the energy minima and obtain the stretch modulus as a coefficient 
of the quadratic term. By giving more weightage to the quadratic term, we 
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have done the higher order fit due to the difficulty in identifying the exact 
Hooke's law region. The term 'sequential' essentially means that the leading 
terms near the energy minima are fit to quadratic term and the rest are fit 
to higher order terms. This procedure is repeated self consistently so that 
the value of 71 is independent of the 7^ for large n. However 71 depends 
very weakly on 72 and 73 (that means 71 has lesser dependence on choice 
of 7„ for n > 1). The energy vs strain curve has a minimum around which 
a quadratic fitting was done and the obtained 71 is the stretch modulus as 
shown in equation O The stretch modulus obtained from a polynomial fit 
to the force-extension curve and also from a sequential fit to the energy- 
strain curve are listed in table [H The values of stretch modulus from the 
above two fitting methods are very similar. From the sequential fit to the 
energy-strain curve, we find that PX-6:5 has a stretch modulus of 1636 pN 
at 246 niM of Na"*" concentration when pulled with 10~^ pN/fs rate, which 
is 30 % more rigid compared to the stretch modulus 1269 pN of B-DNA at 
same concentration and of the same sequence. From the polynomial fit to 
the force-extension curve, the same trend was observed (Table [1]). During 
the pulling the bond lengths and bond angles (Fig. [3]) are changing very 
less whereas the torsion angles (Fig. [3|) are changing almost 100 % with 
respect to the initial values at zero force (See next para for more discus- 
sion). This implies that the backbone helix has large internal dynamics and 
offer less resistance to the applied force. Note that the backbone helix gives 
the structural stability to the DNA molecule. On the other hand, bond 
lengths and bond angles offer more resistance to the applied force. Hence 
the crossover links between helical strands contribute more to the structural 
rigidity of these molecules. Although 30 % of rigidity increase seems small, 
yet this enables one to construct promising DNA nanotechnology devices 
with enhanced rigidity. For the JXi, JX2 topoisomers at the same coun- 
terion concentration we obtain 1515 pN, 1349 pN from sequential fitting 
and 1373 pN, 1521 pN from polynomial fitting respectively. Force-extension 
plot for PX-6:5, JXi and JX2 is shown in Fig. |2(b)| and can see a slight 
change in the slope of linear region. Among all structures, PX-6:5 has large 
stretch modulus both in Na"*" and Mg^+ medium due to the most number 
of crossover points. It was also shown that the optimal design of repeated 
stacks and bundles of nanostructures provide great strength to the molecule 
(js^ ). Experimental reports of stretch modulus for lambda phase DNA is of 
the order of 1000 pN depending on the environmental conditions like ionic 
strength (j24l . l53l ). Though we get the stretch modulus of the same order, 
the helix to ladder transition occurs at a much larger force regime than that 
observed experimentally (j54l . l55l ) due to the higher rates of force employed 
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in our simulation. Experimental results have shown that the ladder trans- 
formation occurs at a much lower (~ 60-70 pN) force than that observed in 
our simulation data. It has been observed that a higher force rate leads to 
higher stiffness of the short DNA (Fig. 2]). This shift of curve and increase 
of slopes in force-extension curves is expected since higher rates means that 
our simulation is far from being reversible, in which case it dissipates energy 
to the environment (since our simulations were constant T = 300 K) and 
higher the irreversibility the higher is the work of dissipation which is given 
by the area under the force-extension curve. A possible reason may also be 
the length of DNA that has been used. Since a 38-mer DNA is too short, its 
contour length being significantly smaller than the persistence length of a 
lambda DNA, its behavior is likely to deviate from that observed for lambda 

DNA Hlsi). 



To gain further microscopic understanding of the structural changes with 
the applied force we have looked at the energetics of the PX/JX DNA as 
a function of the applied force. In Fig. [5] we plot the total internal energy 
as a function of force for PX-6:5, JXi and JX2 DNA structures. The con- 
formational entropy of the structure is dominating for small forces during 
which the change in length is almost negligible. For forces up to 150 pN, 
the energy is decreasing and seems to attain minimum value which corre- 
sponds to a most stable configuration. With further increase in force above 
150 pN, the energy of the DNA molecule increases. This implies that the 
PX/JX DNA structures become thermodynamically unstable under elon- 
gation. Various backbone parameters were calculated when PX-6:5 struc- 
ture is pulled at one end when the structure is neutralized with Na"*^ and 
Mg^"'' counterions. The average bond length (^oa'-p) and average angles 
(0p_O5'_c5' ^^"^ ^C5'-C4'-C3') were changing less than 2 % with applied force 
with respect to the zero force equilibrium values (Fig. [3]). Correspondingly, 
the increase in the energy contribution from bond stretching and angle bend- 
ing is about to 5-10 % compared to the zero force case. Interestingly, all 
the torsion angles (a, /3, 7, 6, e, C) (0) were changing over 100 % (Fig. [3|) 
with respect to the values at zero force, which is causing the structural de- 
formation to great extent. The change in torsion angles is not considerable 
for forces up to 200 pN during which the change in DNA contour length is 
also very less. Dramatic change is observed in the torsion angles when the 
DNA motif elongates to almost twice of its initial contour length at critical 
force. At this stage, except the breaking of few H-bonds (j23l ). we see no 
considerable change in average bond length or average angle. We have used 
geometry measurement based criteria for the H-bond calculation in simula- 
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tion. Generally the H-bond is represented as D— H ■ A where D is the donor 
and A is the acceptor which is bonded to D through the H atom. In the 
case of DNA, D is a N atom and A is either a N or O atoms depending on 
AT and GC base pairing. When the distance between D and A atoms is less 
than 2.7 A and the angle DHA is greater than 130°, we say that the atom A 
is H-bonded to atom D otherwise the H-bond is broken. It is justifiable to 
argue that the backbone atoms in DNA motif were drastically re-oriented 
to give the elastic rod a sudden elongation, a clear signature indicating very 
large change in torsion angles. The torsion energy is also increasing greatly 
compared to all other contributions to the total internal energy of the motif. 
So a closer look at the various energy contributions reveals that, when the 
DNA motif is pulled along the helix axis, there is an increase of about 5-10 
% in the bond stretching and angle bending energy, 15 % increase in van 
der Waals energy but dramatic increase in the torsion energy. Apart from 
all the energy components like bond, angle, van der Waals, Coulomb energy, 
this increase in the torsion energy contributes to rapid increase in the total 
internal energy at a critical force. Instantaneous snapshots of PX-6:5 struc- 
ture at different force in the presence of Na"*" counterions is shown in Fig. [6] 
(see Fig. S4-S13 for instantaneous snapshots of PX-6:5, JXi and JX2 in the 
presence of Na"*" and Mg^"*" at various forces, respectively). As the applied 
force is increased beyond 200 pN, H-bonds in the basepair were broken and 
the DNA overstretches to 1.7 times its initial contour length. 

The original experiments on the PX/JX DNA molecules were performed 
in Mg^"*" buffer (jl8l ). To understand the pulling response of these DNA 
nanostructures with divalent cations we have also done the pulling simula- 
tion of PX-6:5, JXi and JX2 structure in the presence of divalent Mg^"*" ions. 
The counterion concentration in our simulation is close to 160 mM and the 
rate of pulling is 10~^ and 10~^ pN/fs. We have analyzed the force-extension 
spectrum of all DNA structures in Na+ and Mg^^" ions at 10"^ pN/fs (Fig. 
[7]). Also, Fig. S14(a) shows the force-extension behavior of PX-6:5 molecule 
in presence of Mg^"*" ions. In the presence of Mg^+, PX-6:5 have a stretch 
modulus of 1840 pN with a pulling rate of 10~^ pN. In contrast, the stretch 
modulus of B-DNA of same length and sequence as one of the helical do- 
mains of PX-6:5 is 1590 pN. Again we see the enhanced rigidity of the PX-6:5 
motif as was in the presence of Na"*" ions. However, this increase is not as 
dramatic as in the presence of Na^ ions. Possibly the rigidity can be further 
enhanced by increasing the ionic strength. For all the structures we find the 
stretch modulus is more in the presence of Mg^+ than in the presence of 
Na"*" counterions. This is due to the strong phosphate-Mg^"*" coordination 
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that resist the external apphed force. This result is in accordance with the 
experimental results on DNA stretch modulus (jsil ) where it is shown that 
the presence of divalent cations strongly reduce the persistence length of the 
DNA and hence increase the stretch modulus of the DNA. To understand 
the effect of crossover on the rigidity of the DNA motifs in the presence of 
Mg^+ ions we have also done the pulling simulation for JXi and JX2 motifs. 
Fig. S14(b) gives the force-extension behavior of the JXi and JX2 motifs 
and compared with the PX-6:5 in presence of Mg^"*" ions. We get the stretch 
modulus from the sequential fitting of the energy vs strain plot: JXi and 
JX2 have stretch modulus of 1465 pN and 1654 pN respectively. We see 
the same trend with both 10~^ and 10"'^ pN/fs pulling rates as calculated 
from polynomial fitting. To get more of a molecular level picture we have 
plotted the internal energy of the PX/JX structure as a function of pulling 
force as shown in Fig. [5] and we see a similar behavior compared to the en- 
ergy variation observed in the presence of Na"*" ions. Energy is decreasing in 
the small force regime thus making the DNA structures thermodynamically 
more stable. The increase in the energy with further load is rapid compared 
to the case of pulling in Na"*" medium. Increase in the total energy with the 
applied force implies that the PX/JX DNA structure is thermodynamically 
unstable, with some of H-bonds are broken. We give all the stretch modulus 
results in both Na"^ and Mg^+ counterion solution in Fig. S16. For compar- 
ison in Fig. SI and S2 we show the radial distribution functions of the Na''" 
and Mg^+ ions with the OIP and 02P Oxygen of the phosphate backbone. 

3.2 Effect of DNA basepair sequence and rate of pulling 

To investigate the effect of sequence on the structural rigidity of the PX 
crossover DNA motifs we have calculated the force-extension profile for PX- 
6:5 made of only AT or GC basepairs in presence of Na^ counterions. The 
stretch modulus obtained from the sequential fits comes out to be 1592 pN 
for AT sequence and 1780 pN for GC sequence. As expected, a PX struc- 
ture made of only GC sequence is much stiffer than one made of only AT 
sequence due to extra H-bonding in GC basepair. The stretch modulus of 
crossover PX structures made of only AT or only GC is still larger than the 
stretch modulus of 1269 pN of B-DNA double helix of same length with a 
combination of AT and GC sequence. Stretch modulus 71 values for various 
PX-6:5 and JX structures are tabulated in tabled] for comparison. 

It is known from previous theoretical study that the stiffness of the poly- 
mer increase when pulled at faster rate (jisi ) . So when pulling with faster rate 
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the bond strength would increase dynamically and hence stretch modulus is 
expected to increase. Typical AFM pulling rates are of 10^ pN/s whereas the 
slowest pulling rate achieved in our simulation is 10^^ pN/s (i.e., 10~^ pN/fs). 
Therefore due to the higher pulling rates employed in our simulation, the 
obtained stretch modulus is expected to be higher in magnitude compared to 
the results obtained from single molecule experiments. Similarly the magni- 
tude of force where the PX/JX structure elongates roughly twice its initial 
contour length would also be more than that obtained from experiments. 
To see the effect of rate of force on force-extension behavior and stretch 
modulus, we have done simulations at two different pulling rates viz., 10~^ 
pN/fs and 10~^ pN/fs. Fig. 4(a) shows the force-extension curve for PX 
structure with two different pulling rates. It is clear from the plot that the 
force at which the PX structure extends double its initial contour length is 
very high when pulled with 10~^ pN/fs pulling rate compared to 10~^ pN/fs 
pulling rate. This is due to the dynamic stiffening of the H-bonds in PX 
structure when pulled with 10~^ pN/fs compared to 10~^ pN/fs. Breaking 
of H-bonding is the major signature of mechanical deformation of the DNA 
molecule (|57l . l58l ) We have calculated the fraction of surviving H-bonds as 
a function of pulling force at two different pulling rates (Fig. |4(b)[ ). As the 
applied force is increased, the fraction of survived H-bonds decreases and 
at large enough force the fraction of survived H-bonds goes to zero. From 
Fig. |4(b)| it is clear that the fraction of surviving H-bonds goes to zero at 
smaller forces when pulled with 10~^ pN/fs pulling rate compared to 10"'^ 
pN/fs pulling rate. The stretch modulus obtained from energy-strain and 
force-extension curve are listed in table [1] in the presence of Na"*" (Mg^"*") 
counterion medium. 



4 Conclusion 

We have calculated the rigidity of PX DNA molecules in the presence of Na"*" 
and Mg^'^ counterions by directly calculating their force-extension behavior 
under axial stretching. Earlier we demonstrated 
these crossover DNA motifs from the vibrational density of states analysis 
but a quantitative estimate of their structural rigidity was lacking. Now we 
give a quantitative estimate of the stretch modulus of these DNA strictures. 
In the presence of Na"*" ions at a counterion concentration 246 mM the stretch 
modulus of the PX-6:5 structure is almost 30 % more than that of normal 
B-DNA double helix of same length and having sequence of one of the double 
helical domain of PX-6:5. The computational cost of these calculations is 



20, M, m) the rigidity of 
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enormous, thus restricting us to use intermediate pulling rates of 10~^ and 
10"'^ pN/fs. To understand the effect of crossovers on the stretch modulus 
of these DNA motifs we have also calculated the force-extension profile of 
the JX1/JX2 motifs in the presence of Na^ counterions. We find that JXi 
has stretch modulus of 1515 pN which is slightly smaller than the stretch 
modulus of PX-6:5 (1635 pN). JX2 has a stretch modulus of 1349 pN which 
is 286 pN smaller than the stretch modulus of PX-6:5. When the DNA is 
pulled, among all contributions to the total energy of the DNA, there is a 
dramatic increase in torsion energy. This increase in torsion energy is due to 
a very large change in different torsion angles (Fig. [3]) which cause the DNA 
to destabilize with increased force. Interestingly there is almost no change in 
various bond distances or angles as a function of force (Fig. [3]). The similar 
behavior was observed in the presence of divalent Mg^"*" ions. In the presence 
of Mg2+ we find the stretch modulus of PX-6:5, JXi and JX2 to be 1875 pN, 
1465 pN and 1654 pN, respectively. PX-6:5 has the highest stretch modulus 
than any other DNA motif as in the case of Na^ medium. Interestingly all 
structures in Mg^"*" medium have more stretch modulus compared to Na"*" 
medium due to large electrostatic screening arising from the divalency of 
Mg^"*" ions. This could be due to the fact that presence of Mg^+ ions gives 
extra stability to the structure, owing to the strong coordination of Mg^"*" 
with the phosphate atoms of the two double helical domains. We have also 
studied the effect of the rate of force on the rigidity of DNA motif and found 
that the increased rate of force enhances the rigidity of the structure. 

5 Acknowledgements 

We acknowledge Department of Science and Technology (DST), Government 
of India for financial support. PKM also thanks Alexander von Humboldt 
foundation for sponsoring his visit to Technical University Munich where 
part of the work was done and Roland Netz for help with the sequential 
energy fitting and valuable comments. We are also grateful to Prof. Ned 
Seeman for a critical reading of the manuscript and valuable suggestions. 
MS thanks University Grants Commission (UGC), India for senior research 
fellowship. 



Structural Rigidity of PX and JX DNA 



15 



References 

1. Seeman, N., 2003. DNA in a material world. Nature 421:427-431. 

2. Seeman, N., 2001. DNA nicks and nodes and nanotechnology. Nano 
Lett. 1:22-26. 

3. Cohen, S., A. Chang, H. Boyer, and R. Helling, 1973. Construction of 
biologically fmictional bacterial plasmid in-vitro. Proa. Natl. Acad. Sci. 
USA. 70:3240-3244. 

4. Seeman, N., 1982. Nucleic- Acid junctions and lattices. J. Theor. Biol. 
99:237-247. 

5. Chen, J., and N. Seeman, 1991. Synthesis from DNA of molecule with 
the connectivity of a cube. Nature 350:631-633. 

6. Zhang, Y., and N. Seeman, 1994. Construction of a DNA-truncated 
octahedron. J. Am. Chem. Soc. 116:1661-1669. 

7. Du, S., B. Stollar, and N. Seeman, 1995. A synthetic DNA molecule in 
3 knotted topologies. J. Am. Chem. Soc. 117:1194-1200. 

8. Du, S., and N. Seeman, 1994. The construction of a trefoil knot from a 
DNA branched junction motif. Biopolymers 34:31-37. 

9. Mao, C, W. Sun, and N. Seeman, 1997. Assembly of Borromean rings 
from DNA. Nature 386:137-138. 

10. Yurke, B., A. Turberfield, A. Mills, F. Simmel, and J. Neumann, 2000. 
A DNA-fuelled molecular machine made of DNA. Nature 406:605-608. 

11. Buehler, M. J., and T. Ackbarow, 2007. Fracture mechanics of protein 
materials. Materials Today 10:46-58. 

12. Seeman, N., H. Wang, X. Yang, F. Liu, C. Mao, W. Sun, L. Wenzler, 
Z. Shen, R. Sha, H. Yan, M. Wong, P. Sa-Ardyen, B. Liu, H. Qiu, X. Li, 
J. Qi, S. Du, Y. Zhang, J. Mueller, T. Fu, Y. Wang, and J. Chen, 1998. 
New motifs in DNA nanotechnology. Nanotechnol. 9:257-273. 

13. Fu, T., and N. Seeman, 1993. DNA double-crossover molecules. Bio- 
chemistry 32:3211-3220. 

14. Sa-Ardyen, P., A. Vologodskh, and N. Seeman, 2003. The flexibility of 
DNA double crossover molecules. Biophys. J. 84:3829-3837. 



Structural Rigidity of PX and JX DNA 



16 



15. Winfree, E., F. Liu, L. Wenzler, and N. Seeman, 1998. Design and 
self-assembly of two-dimensional DNA crystals. Nature 394:539-544. 

16. Holliday, R., 1964. Mechanism for gene conversion in fungi. Genet. Res. 
5:282-&. 

17. Li, X., X. Yang, J. Qi, and N. Seeman, 1996. Antiparallel DNA double 
crossover molecules as components for nanoconstruction. J. Am. Chem. 
Soc. 118:6131-6140. 

18. Shen, Z., H. Yan, T. Wang, and N. Seeman, 2004. Paranemic crossover 
DNA: A generalized Holliday structure with applications in nanotech- 
nology. J. Am. Chem. Soc. 126:1666-1674. 

19. Yan, H., X. Zhang, Z. Shen, and N. Seeman, 2002. A robust DNA 
mechanical device controlled by hybridization topology. Nature 415:62- 
65. 

20. Maiti, P., T. Pascal, N. Vaidehi, and W. Goddard, 2004. The stability 
of Seeman JX DNA topoisomers of paranemic crossover (PX) molecules 
as a function of crossover number. Nucleic Acids Res. 32:6047-6056. 

21. Maiti, P., T. Pascal, N. Vaidehi, J. Heo, and W. Goddard, 2006. 
Atomic-level simulations of Seeman DNA nanostructures: The parane- 
mic crossover in salt solution. Biophys. J. 90:1463-1479. 

22. Spink, C. H., L. Ding, Q. Yang, R. D. Sheardy, and N. C. Seeman, 
2009. Thermodynamics of Forming a Parallel DNA Crossover. Biophys. 
J. 97:528-538. 

23. Santosh, M., and P. K. Maiti, 2009. Force induced DNA melting. J. 
Phys.: Condens. Matter 21:034113. 

24. Bustamante, C., J. Marko, E. Siggia, and S. Smith, 1994. Entropic 
elasticity of lambda-phase DNA. Science 265:1599-1600. 

25. Smith, P., and B. Pettitt, 1996. Ewald artifacts in liquid state molecular 
dynamics simulations. J. Chem. Phys. 105:4289-4293. 

26. Allemand, J., D. Bensimon, R. Lavery, and V. Croquette, 1998. 
Stretched and overwound DNA forms a Pauling-like structure with ex- 
posed bases. Proc. Natl. Acad. Sci. USA. 95:14152-14157. 



Structural Rigidity of PX and JX DNA 



17 



27. Rief, M., H. Clausen-Schaumann, and H. Gaub, 1999. Sequence- 
dependent mechanics of single DNA molecules. Nat. Struct. Biol. 6:346- 
349. 

28. Vologodskii, A., 1994. DNA extension under the action of an external 
force. Macromolecules 27:5623-5625. 

29. Marko, J., and E. Siggia, 1995. Stretching DNA. Macromolecules 
28:8759-8770. 

30. Smith, S., Y. Cui, and C. Bustamante, 1996. Overstretching B-DNA: 
The elastic response of individual double-stranded and single-stranded 
DNA molecules. Science 271:795-799. 

31. Cluzel, P., A. Lebrun, C. Heller, R. Lavery, J. Viovy, D. Chatenay, and 
F. Caron, 1996. DNA: An extensible molecule. Science 271:792-794. 

32. Fu, H., H. Chen, J. F. Marko, and J. Yan, 2010. Two distinct over- 
stretched DNA states. Nucleic Acids Res. 38:5594-5600. 

33. Fu, H., H. Chen, X. Zhang, Y. Qu, J. F. Marko, and J. Yan, 2011. Transi- 
tion dynamics and selection of the distinct S-DNA and strand unpeeling 
modes of double helix overstretching. Nucleic Acids Res. 39:3473-3481. 

34. Chen, Hu and Yan, Jie , 2008. Effects of kink and flexible hinge defects 
on mechanical responses of short double-stranded DNA molecules. Phys. 
Rev. E 77:041907. 

35. Rouzina, I., and V. Bloomfield, 2001. Force-induced melting of the DNA 
double helix. 2. Effect of solution conditions. Biophys. J. 80:894-900. 

36. van Mameren, J., P. Gross, G. Farge, P. Hooijman, M. Modesti, 
M. Falkenberg, G. J. L. Wuite, and E. J. G. Peterman, 2009. Unraveling 
the structure of DNA during overstretching by using multicolor, single- 
molecule fluorescence imaging. Proc. Natl. Acad. Sci. USA. 106:18231- 
18236. 

37. Cocco, S., J. Yan, J. Leger, D. Chatenay, and J. Marko, 2004. Over- 
stretching and force-driven strand separation of double-helix DNA. 
Phys. Rev. E 70:011910. 



38. Lai, P., and Z. Zhou, 2003. B- to S-form transition in double-stranded 
DNA with basepair interactions. Physica A 321:170-180. 



Structural Rigidity of PX and JX DNA 



18 



39. Tung, C, and E. Carter, 1994. Nucleic-acid modeling tool (NAMOT) - 
an interactive graphic tool for modeling nucleic-acid structures. Comput. 
Appl. Biosci. 10:427-433. 

40. Case, D. A., T. A. Darden, T. E. Cheatham, I., C. L. Simmerling, 
J. Wang, R. E. Duke, R. Luo, K. M. Merz, B. Wang, D. A. Pearlman, 
M. Crowley, S. Brozell, V. Tsui, H. Gohlke, J. Mongan, V. Hornak, 
G. Cui, P. Beroza, C. Schafmeister, J. W. Caldwell, W. S. Ross, and 
P. A. Kollman, 2004. AMBER8. 

41. Duan, Y., C. Wu, S. Chowdhury, M. Lee, G. Xiong, W. Zhang, R. Yang, 
P. Cieplak, R. Luo, T. Lee, J. Caldwell, J. Wang, and P. Kollman, 
2003. A point-charge force field for molecular mechanics simulations of 
proteins based on condensed-phase quantum mechanical calculations. J. 
Comput. Chem. 24:1999-2012. 

42. Wang, J., P. Cieplak, and P. Kollman, 2000. How well does a restrained 
electrostatic potential (RESP) model perform in calculating conforma- 
tional energies of organic and biological molecules? J. Comput. Chem. 
21:1049-1074. 

43. Jorgensen, W., J. Changrasekhar, J. Madura, R. Impey, and M. Klein, 
1983. Comparison of simple potential function for simulating liquid 
water. J. Chem. Phys. 79:926-935. 

44. Aqvist, J., 1990. Ion water interaction potentials derived from free- 
energy perturbation simulations. J. Phys. Chem. 94:8021-8024. 

45. Baker, N., D. Sept, S. Joseph, M. Hoist, and J. McCammon, 2001. 
Electrostatics of nanosystems: Application to microtubules and the ri- 
bosome. Proc. Natl. Acad. Sci. USA. 98:10037-10041. 

46. Baker, N., 2005. Improving implicit solvent simulations: a Poisson- 
centric view. Current Opinion in Structural Biology 15:137-143. 

47. Ryckaert, J., G. Ciccotti, and H. Berendsen, 1977. Numerical- 
integration of cartesian equations of motion of a system with constraints 
- molecular-dynamics of N-alkanes. J. Comput. Phys. 23:327-341. 

48. Evans, E., and K. Ritchie, 1997. Dynamic strength of molecular adhe- 
sion bonds. Biophys. J. 72:1541-1555. 



Structural Rigidity of PX and JX DNA 



19 



49. Hugel, T., M. Rief, M. Seitz, H. Gaub, and R. Netz, 2005. Highly 
stretched single polymers: Atomic-force-microscope experiments versus 
ab-initio theory. Phys. Rev. Lett. 94:048301. 

50. Seol, Y., J. Li, R C. Nelson, T. T. Perkins, and M. D. Betterton, 2007. 
Elasticity of short DNA molecules: Theory and experiment for contour 
lengths of 0.6-7 mu m. Biophys. J. 93:4360-4373. 

51. Wiggins, P. A., T. Van der Heijden, F. Moreno-Herrero, A. Spakowitz, 
R. PhiUips, J. Widom, C. Dekker, and P. C. Nelson, 2006. High flexibil- 
ity of DNA on short length scales probed by atomic force microscopy. 
Nature Nanotechnol. 1:137-141. 

52. Ackbarow, T., and M. J. Buehler, 2009. Alpha-helical protein domains 
unify strength and robustness through hierarchical nanostructures. Nan- 
otechnology 20:075103. 

53. Baumann, C, S. Smith, V. Bloomfield, and C. Bustamante, 1997. Ionic 
effects on the elasticity of single DNA molecules. Proc. Natl. Acad. Sci. 
USA. 94:6185-6190. 

54. Konrad, M., and J. Bolonick, 1996. Molecular dynamics simulation of 
DNA stretching is consistent with the tension observed for extension 
and strand separation and predicts a novel ladder structure. J. Am. 
Chem. Soc. 118:10989-10994. 

55. Lebrun, A., and R. Lavery, 1996. Modelling extreme stretching of DNA. 
Nucleic Acids Res. 24:2260-2267. 

56. Saenger, W., 1984. Principles of Nucleic Acid Structure. Springer Ver- 
lag, New York. 

57. Keten, S., and M. J. Buehler, 2008. Asymptotic Strength Limit of 
Hydrogen-Bond Assemblies in Proteins at Vanishing Pulling Rates. 
Phys. Rev. Lett. 100:198301. 

58. Keten, S., and M. J. Buehler, 2008. Geometric confinement governs the 
rupture strength of H-bond assemblies at a critical length scale. Nano 
Lett. 8:743-748. 

59. Maiti, P. K., T. A. Pascal, N. Vaidehi, and W. A. Goddard, HI, 2007. 
Understanding DNA based nanostructures. J. Nanosci. Nanotechnol. 
7:1712-1720. 



Structural Rigidity of PX and JX DNA 



20 



Table 1: Stretch modulus 71 (pN) from sequential fit to the energy-strain 
curve and polynomial fit to the force-extension curve in the presence of Na"^ 
(Mg^+) counterions with different rates offeree of 10~^ and 10"^ pN/fs (Fig. 
S15). Note that the values in brackets corresponds to Mg^~^ case. 



71 for Na+ (Mg^+) (pN) 
Sequential Fit Polynomial Fit 



DNA 


c (mM) 


10"^ 


10-^ 


10-* 


10-^ 






(pN/fs) 


(pN/fs) 


(pN/fs) 


(pN/fs) 


B-DNA 


254(140) 


1269 


1358 


1578(1591) 


1644 


PX-6:5 


246(160) 


1636(1875) 


1737 


1772(1841) 


2322 


JXi 


246(160) 


1515(1465) 


1558 


1374(1404) 


2176 


JX2 


246(160) 


1349(1654) 


1480 


1521 


2051 


AT65 


246 


1592 


1731 


1378 


1981 


GC65 


246 


1780 


1862 


1546 


2138 
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Figure Legends 

Fig. m 

Atomic level structure of PX/JX DNA molecules: (a) The basepair se- 
quences used in the generations of PX-6:5, JXi and JX2 crossover molecules. 
Initial structure of (b) PX-6:5, (c) JXi and (d) JX2 used in our pulling sim- 
ulation. 

Fig. El 

Force-extension curves: (a) Force-extension behavior of PX-6:5 and B-DNA 
molecule in the presence of Na"*" counterions. (b) Force-extension behavior 
of PX-6:5, JXi and JX2 in presence of Na^ ions. The rate of force is 10"^ 
pN/fs. 

Fig. m 

Variation of backbone parameters with force: 03'-P bond distance, P-05'- 
C5' angle and C5'-C4'-C3' angle as a function of force applied on PX-6:5 
in Na"*" (black) and Mg^+ (blue dashed line) medium. There is very little 
change in the bond distance and the angles with the increase in force. Vari- 
ation of torsion angles alpha (a), beta (/3), gamma (7), delta (6), epsilon 
(e) and zeta (C) as a function of force applied at the end of PX-6:5 in Na"*" 
(black) and Mg^"*" (blue dashed line) medium. There is a very large change 
over 100 % of the equilibrium zero force limits (dotted lines) in all the torsion 
angles with the force applied. 

Fig. m 

Effect of rate of pulling on PX-6:5 structure: (a) Force-extension curve at 
two different pulling rates. Higher the pulling rate steeper is the response 
and hence large stretch modulus, (b) Fraction of survived H-bonds as a 
function of pulling force at different pulling rate. At high pulling rate melt- 
ing/breaking of H-bonds occurs at larger magnitude of pulling force. 

Fig. m 

Energy of various PX/JX structures as a function of pulling force in Na^" 
and Mg^~^ counterions. The rate of pulling is 10~^ pN/fs. The minimum in 
the energy curve corresponds to the most stable structure during pulling. 
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Fig. m 

Instantaneous snapshots of the PX-6:5 structure at various forces (a) 200 
pN, (b) 400 pN, (c) 600 pN, (d) 800 pN and (e) 1000 pN in the presence 
of Na+ counterions. As the applied force is increased beyond 200 pN, H- 
bonds in the basepair were broken where the extension in the DNA is almost 
double its initial contour length. 

Fig. m 

Force-extension spectrum analysis. Upper color bar indicates the strain in- 
crease percentage on applying the external force on the various DNA struc- 
tures. The extension in Mg^"*" medium is less than the extension in Na"*" 
medium which implies the enhanced rigidity of the PX/JX DNA molecules 
in Mg^"*" medium. 
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